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Initial proposed themes 
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‘ The proposal, from the french Meteorologie Nationale, 
for participation in SKYLAB EREP Experiment was formally adopted 
by NASA, in October 1972 (Letter SR 008, dated October 6, 1972) , 


The proposal included three themes 
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The theme N° 1 followed the lines of a current scien- 
tific project, lead by french meteorologists, investigating 
luesoscale cloud features physical mechanisms, and using photo- 
graphs taken by constant level high altitude balloons, simulta- 
neously vrith thermodynamic and kinematic informations given by 
an associated ground network. 




Its main goal was to improve a model describing the 
physical relationship between the convective cloud rows and the 
wind profiles, and to check its rehability. 


The theme N° 2 used the former demonstration that micro- 
wave attenuation and diffusion through clouds could give infor- 
mation ::ibout the water content in cloud depth. It had seemed 
possible to infer, from such radiometric data, the occurence of 
precipitation under clouds over the ocean, EREP data would give 
an opportunity to check the mathematical model rehability. 


The theme N® 3 wanted to determine the snow cover, its 
temperature, the beginning of the melting and to research sta- 
tistical relations between superficial characteristics of the 
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snow and its spectral signature. 

Proposed associated experiments 
Theme_N®_l 

There would have been a joint experiment conceived as 
follows : EREP would give a detailed picture of cloudiness, both 
by its photographic facility S 190, and its IR sensor S 192, ope- 
rating in the 10 to 12,5 micrometer atmospheric windov/, over 
south-west of FRANCE. Simultaneously, according to satellite ephe- 
meris, a very dense ground network of radiowind, radiosonde and 
radar stations, specially settled in the same region, and supple- 
mented by three specially equipped aircrafts, would operate and 
give information on the thermodynamic and kinematic structure of 
the observed air mass (fig. 1,2,3). 

Theme_N®_2 ' 

The experiment would make use of EREP data r CCT, at va- 
rious times of the same part of the Atlantic Ocean, for various 
meteorological situations, without cloud, with thin or thick 
precj-pitation. A special mode of operation, for S 193, was asked 
and accepted. The ocean temperature v/ould be simultaneously che- 
cked by NIMBUS HRIR measurements. The cloud droplets distribution 
could be obtained by in situ mesurements; twenty-five flight hours 
of a french Navy plane "NEPTUNE" were scheduled for this purpose 
(fig. 4) . 

^ : ■ ' h. 

^ Ground mesurements of snow state and temperature were t ■ 

foreseen, simultaneously with SKYLAB S 192 data (fig. 5). *; 

' .'t 

How the experiments worked ■ ^ 

■ ^ — — ■ I ' 

Theme_N®_l .i 

' 

The chbosen date for the experiment was September 17, 1973; 
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A very important ground and air network v/as in position, inclu- 
ding 12 complete radio sounding stations, 2 radar stations, l 
plane used for upper air measurements and 2 high altitude MIRAGE 
IV planes were used for taking pictures of nebulosity along one 
choosen orbit path. The track displacement of 2®46' East made 
this netv/ork completely useless. 

Theme N® 2 

mm mm mm mtm 'mm mm 

The experiment, which was scheduled for September 1973 

i * 

(SL 3) , had to be cancelled, because of an ill-functioning of the 
S 193 SKYLAB apparatus. 

Theme_N®_3 

The experiment, which was scheduled for the end of May 
1973, had to be cancelled, because of the slippage of the SKYLAB 
II launching and the actual fact that the snow melting is to be 
completed by raid- June over the french ALPES. 

Final results 

As we just said, there were no result for experiment 
505802 and 505804, which had to be cancelled, and for the first 
experiment, the experiment 505801, there was also no. result for 
the choosen date of September 17. 

^ But, as it is mentionned above, the main goal of the 

theme N° 1 was to check and to improve amodel describing the 
physical relationship between the convective cloud rows and some 
thermodynamic parameters. Without a special network, it has been 
possible, for this purposey to work' .on a very interesting cloud 
streets pictures, taken over FRANCE, by SKYLAB S 190, the 19th 
of September 1973, 

The result of this work, a document about "The possibili- 
ty of evaluating wind profiles from satellite data" was sent to 
you March 24, 1975. I should appreciate if we could receive some 
kind of appreciation or criticism about this paper. 
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PARTICIPATION IN THE SKYLAB PROGRAM 


The possibility of evaluating vertical wind profiles from satellite data 

wtiLL e r<. 

ABSTRACT 

In the presence of convective movements, a method is proposed for evalua' 
ting the vertical variation of the horizontal wind. The profile is 
determined from the following data; on the one hand, the wind and 
the air temperature at the upper and lower boundaries of the convective 
layer; on the other hand, the value of a so-called characteristic 
dimension enabling the determination of thtee wavelengths which define 
the general morphology of the convective cell. 

The author, after discussing the establishment of the vertical wind 
equation, proposes a sinusoidal solution for it, which is appropriate • 
to the phenomenon investigated. This solution requires a clearly 
defined profile of the horizontal wind. The morphology of the cell 
is derived from its characteristic dimension by a minimization of 
the internal dissipation of friction energy, the calculation being 
made for six types of cell. The author provides details concerning 
the numeric analytical methods employed and supplies a general 
calculation flowchart- Finally, tables are presented covering widely 
varying regions, corresponding either to cloudstreets or to cells, 
depending on the vertical thermal gradient and the value of the 
characteristic dimension. 
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This method vas developed during an in-flight experiment on Skylab III, 
which permitted accurate measurements if a cloudstreet observed over 
France. A comparison was made of the profile calculated by this 
method and that measured by the nearest wind sounding. 

As a consequence of the initial satisfactory results thus obtained, 
the author proposes to systematically check this method on a 
sufficient number of radiosoundings, the thickness of the convective 
layer: being taken as the characteristic dimension of the convective 
cell. Subsequently, should the results prove satisfactory, the method 
may be applied to satellite pictures of cloudstreets, taking the 
inter-street spacing as the characteristic dimension, and calculating 
the thermal gradient as a function of the service temperatures 
supplied by the V.H.R.R. 





VERTICAL WIND EQUATION 

Throughout this study, the symbols employed by J.P. Kuettnor 
in his article published in 1971 (1) will be employed as far 
as possible. 


It is worthwhile starting with an initial reminder of the 
establishment of the vertical wind equation. In the initial 
stage, the turbulent viscosity and conductivity coefficients 
are ignored. The basic equations, reflecting a situation of 
stationary state convection, can be written as follows: 
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with g “ I ^ I ' 

Although this model of equation is generally acceptable for 
this type of problem, it is worthwhile elucidating some of 
its underlying hypotheses. 
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,No account is taken of the Coriolis force, which is negligible 
at the scale of circulation of convection with respect to the 
local horizontal pressure gradient. The continuity equation 
is represented by non~divergence of the wind vector. This 
implies an acceptance of Boussinesq's so-called incompressi- 
bility hypothesis. This approximation postulates that at a 
given level, the density of the air varies only with tempera- 

t * 

ture. 

By giving, V A, and to the symbols A»A,A',A*', the same relative 
meanings as Kuettner (1), A representing the perturbation term 
and A' the first derivative with respect to z, we have the 
following equivalents: 


5 (5 — — « — — (5 - 

= W = -gj (u,v,p,p,T) = -gp- (u,v,p,p,T) = 0 


Hence: 

A = - 
§ - 6z 

and, in accordance with Boussinesq's approximation, 

± - 1 
. Y T 

Inlthese conditions, by employing the operator h, such that 


and by linearizing, namely, ignoring the minor terms of the 


type 
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the initial system can be written; 
1 


h (u) + u’w - 
Ih (v) + v'w = 


h (w) + J T = 
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+ w’ = 0 
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h (T) + (7^ + T')w = 0 

SL 


2.5 


Starting with these basic equations, this is the procedure to 
be followed in order to obtain the vertical wind equation. The 
operators and g - - are applied to 2,1 and 2.2 respectively; 

hence, after linearization; 
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and, by addition, by employing 2.4, being the Laplacean operator 
and , the horizontal Laplacean operator 


- » |l ^ s? =■ - ='^5 p’ 


The substitute operator Vj is applied to 2.3 with h 


(4)' 


(h(w)) - § Vj (T) = -fVj (p’) 

T P 

I 

and, by subtracting (5) from (4), 

- Vt (h(w)) + 5 V? (T) + + v’* = 0 
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(5) 
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Finally, by applying the operator h to (6) while elucidating 
h (T) in accordance with 2.5 

-h^ (V2(w)) - I (y^ + T’) Vf (w) + h (u" + v" 3 ^) =0 (7) 

This is the vertical wind equation which, in this form, is an 
equation of the fourth order. 


SOLUTION OF THE VERTICAL WIND EQUATION IN THE PRESENCE OF 
CONVE.CTIVE CELLS 

The convective cells can be considered to be a stationary phenomenon 
spatially repetitive, comprising at the interior of each cell a 
sinusoidal vertical wind variation, along the vertical axis and 
two horizontal axes perpendicular to each other. The vertical 
wind field hence admits three wavelengths L^, Ly and L^, enabling 
the introduction of the numbers: 


2tt 
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As an expression of this field, we 


w = Wjj cos (lx) cos (my) sin (nz) 


( 8 ) 


0SI6INAL PAGE IS 
D1 POOR qUAUW 


In the presence of clouds, the maximum ascendancy w^^ occurs 
immediately below the clouds and the maximum descendancy -w^ 
occurs between the clouds. In the presence of cloudstrcets, 

Ljj is the spacing between the centers of two Consecutive 
clouds in a single street, Ly is twice the spacing between the 
axes of two streets, and L^ is twice the thickness of the 
convective layer. The following six types of cell were studied: 
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Lx = Lz 

2Lx = Lz 

2Lx - 2Lz ■ 

4Lx = 2Lz 

l/i Lx = 2/2 Lz 
Lx * Lz 


These provide for a gradual passage from aerological scale to. 
the mesoscale, by a gradual increase of horizontal dimensions 
with respect to the vertical dimension. 


Types 2, 3 and 4 may correspond to cloudstreets , and types 5 
and 6 to open and closed scales of the mesoscale. The structure 
suggested for the vertical wind field appears to be most 
reasonable, particularly in the presence of cloudstreets. It is 
different from that adopted by J.P. Kuettner in his 1971 article. 


By introducing structure (8) in the vertical wind equation (7) 
it will be shown that this may be a solution, provided that one 

. i ' ' . ’ ' . 

Accepts tWATt' cons train ts the parameters of this equation 

i ' ■ - . 

These constraints constitute the interesting part of the model, 

by conditioning the vertical variation of the horizontal wind. 

We are thus faced with a reverse method for determining the. 

♦ 

■ ■ ■ ■ ■ ■ . 

wind profile. 
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PRINCIPLE OF THE DETERMINATION OF THE WIND PROFILE 
Assuming the following auxiliary variables: 

,2f 2- 2 

d=l +m.+n 
X = sin (lx), X* = cos (Ix) 

Y = sin (ny) , Y’ - cos (my) 

Z - sin (nz)', Z* = cos (nz) • . 

If ve replace W, the Laplacean operators, the operator h and 
. the partial derivatives of w by their expressions as a function 
of X, Y, Z, X', Y', Z', we obtain: 

d^h^(w) + ;§ (y^ + T)(l^ + m^)w + h(u" ^ + v" ^) =0 (9) 

with 

d^h^(w) - d^W^j |j“ l^<j? - v^)X’Y’Z + 21 m u v X Y z] 

- (Y^ + Y')(l^ + m^)w = ^ (y + T')(l^ + m^) X'Y'Z 

f ® ^ M a 

h(V." ■“ + v" -y) - V^j |](~ l^u u - m^v V ”)X’Y*Z+lra(v u *'+u V ")XY^ (10) 

The vertical wind equation can then be written 

lm(2d^u V + u v" + v u")XY + [^“(Y +T’)(l^+m^) - 1^^ u'*-m^^ V - 

/(l^u^ + m^v^r|x'Y' = 0 

; ; ■ - / 

Sihce this equation must be satisfied at all points x,y, it must 
also be satisfied for each of the two pairs XY and X'Y'. We must 
thys have a vertical wind structure such that at any level z, the 
two horizontal wind components u and v satisfy the system: 

, lU:., 
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■2 u V + u v" + V u" = 0 

(Y^ + T )(1^ + m^) ” u u" ~ V v" - d^(l^u^ + m^v^) =0 (M) 

T ^ 


In other words, by writing the vertical thermal gradient with 
its usual notation Y* nnd taking as the air temperature at 
the ground: 


T7» = f g (1^ (Y - Y o-L , _ ^2” „ 

j_1(f - \z)(- l^u^ + m2v2 J 


g(l^ •»• m^)(T - Ygu ) ^ - d^ 

(T^ - z) (l^u'^ - m2v2) 


02 ) 


Note that for an adiabatic gradient, we have 


—II 

V 


•<«v % ^ tiWf,. 

In other words, a sinusoidal wind variation.vistdfcja.-Jimaxffesium 

tht Ativ kA«itu •{- W 

<fi^':fcltE>niiddli?'<<rf<>tho-.-eai^veG^4ye6>^ 


Within a convective cell, system (12) can be accepted as a 
necessary constraint governing the wind profile. It can be 
noted that in the case of streets which are indefinitely 
oriented in the wind direction, namely, in the case of two- 
dimensional rolls, system (11) is reduced to 


-4- (Yjt + T')m^ =0 

^ A. 1 

leading to the need for an adiabatic thermal gradient. 


05) 


in 
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NUMERICAL CALCULATION OF THE WIND PROFILE 

To process numerically equation (12), we approximate the 
secondary derivatives u" and v" by finite differences of 
the type: 


uV = “i^l ~ ^ Vl 

^ . Az^ 


v.^, - 2v. + V. , 

V” s: ^ l~l 

i a_2 


Az' 


where u and v are the wind components at level iAz. As a 
value of Az, it is sufficient to take a number of the order 
of Lz/ 1000. In practice, one generally takes Az = 10 m. 

By writing: 

X _ 

we obtain the recurrence formulas 

U. _ A;., A> _ «: 

A - ^ X. ^ V ay A. -I • 


A;..- Ay 


A - 'I 


- 1 t t. t 

I AA. ^ - yy\ '\r 






A - A 


( 16 ) 


T^he two horizontal wind components u^ and v^ at level iAz, along 
the two perpendicular direction L and L , may thus be calculated 

_ * y 

by iteration, as a function of U£_j» '^i-i* '^i-l* 


II. 


The wind can thus be calculated at all levels, if it is known 
at the two first levels, level 1, the ground level, i.,nd level 2, 
the Az level. 

However, it is easier to know the wind at level I 2 / 2 , the top 
of the convective layer, than at level 2. The wind at level 2 
may be obtained by a rapidly convergent calculation of succes- 
sive approximations, based on a linear error interpolation, of 
which the algorithm is given below. 

The procedure suggested for obtaining the wind profile, based 

on equations (16), is only consistent if the denominator 

,22 22 , ... . . . 

1 u^ + m v^ always maintains its original sign, which must 

be checked at each step. The wind is only calculated inside 

the convective layer '^/ 2 \ it is printed at 100 meter intervals. 


5. ALGORITHM FOR CALCULATING THE WIND AT LEVEL A z 


Let the values of the wind components at the ground be URj , VRj 
and at the top UR^j and VRjj. For level 2, corresponding to 
level Az, we first assume two arbitrary values, for example t. 


UG, = UR, + A X 10 
1 I z 


VGj = VRj + Az X 10 


The iterative calculation from URj , UGj, VRj and VGj leads to 
level H corresponding to at values UCj and VCj. We thus 

■ 
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know the errors 
P| = URj - UCj 

Qj = VRj - VC, 

These cannot be acceptable unless 

P| + Q| < 0.01 (17) 

t * ' 

If not, we can calculate the successive terras of a guess field 
of the values of the wind components at level 2, UGjj and VG^^, 
to finally arrive, by employing the iteration equations (16), 
at the wind component values at level H satisfying equation (17). 

In order to do this, we take the preceding terms as the succes- 
sive values of the guess field, corrected by linearly extrapolating 
the simultaneous influence of a very slight variation in UG and VG 
on the two components of the error at the top and Q^j. 

■ I 

*’n “ ■ '’’'N 

F6r example, VG,^ and (UG^, + 0.0001) produce differences dp, in 
and dqj in and (VG,^ + 0.0001) and UG,^ produce differences 
dp 2 in Pjj and dq 2 in Qjj. 

We have the equations: | y 


X dp, + y dp2 = P,^ 
X dq, + y dq 2 = 
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From this we deduce: 

X = ^ ^ ■^‘Ii ^U- 

<^r,i ‘‘ii " I'*- K 

and the iteration formulas 
UG(N). = UG(fJ-I) + 0.001 X . 

VG(N) = VG(N“1) +0.001 y (18) 

We leave the algorithm for one of the two following conditions, 
stated arbitrarily: 

P(N)^ + Q(N)^ <0.01 
N = 10 

In fact, this system is very rapidly convergent and one generally 
obtains the wind at level Az for N *= 2 or 3. 


DETERMINATION OF THE TYPE OF CELL AND OF THE CORRESPONDING 
WIND PROFILE 

Up to this point, we have ignored the friction forces. They will 
now be accounted for by assuming that the type of cell which 

■ : i . ■ 

should be retained from among the six types considered will be 
that corresponding to a minimization of these forces. 


u. 


.If we assume that the friction forces have the form V 


15 

6z^ 


with a turbulent kinematic viscosity coefficient V which is 
constant, the friction force per unit weight is proportional 
at any level to V . 


For each type of cell, the following sum is calculated 

‘I Vi -■ v; ■ 


(19) 


i = 10,20,30,. .. ,L /2-10- 

z 


This .represents, at the thickness of the convective layer, a 
quantity proportional to the energy dissipated per unit time 
by the friction forces. It should be noted that, owing to 
the fact that each term is a scalar product of two vectors of 
which the components are known, its calculation is especially 
easy to program. Finally, the type of cell retained is that 
corresponding to the minimum value of this sum. 


Me thus obtain the wind profile at the interior of the convective 
layer, for a given characteristic dimension, whether this is the 
thickness of the convective layer, derived from a radio sounding, 
or the inter-street spacing, measured from a satellite picture, 
and from the information concerning the movement and temperature 
of the air at the ground level and at the top of the convective 
layer. J 


Th|e following page shows the table giving the type of cell to be 
re'tained, and consequently enabling the calculation of the corres- 
ponding wind profile, as a function of y and L^, y varying from 

-A _4 „i 

50.10 to 100,10 deg.m and L varying from 2 to 100 km, for 

Pr n 




TYPES OF CELL AS A FUNCTION OF 


L y AND Y 


0 - No solution 

1 - L2 - Ly L^ - , 

2 “ L^ ” Ly Ljj^ “ Ly/2 

3 - = Ly/2 Lj^ = Ly/2 

4 - L^ . Ly/2 Lx - Ly/4‘ 

5 - = Ly/2/f Lx=Ly/2v^ 

6 - L^ = Ly/2*^ Lx=l^/4yf 


100000 vvt 
90000 
80000 
70000 
60000 
50000 
40000 
30000 
28000 
26000 
24000 ! 
22000 
20000 
18000 
16000 
14000 
12000 
10000 
9000 
8000 
7000 

60d0 

5060 

4000 

3000 

2000 


70 75 


lOO. 10 ^ deg, in * 
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a ground wind at 270® and 8 ms ^ , a wind at the top of the 
convective layer at 280® and 9 ms \ a street direction of 
280, corresponding to an x axis of 250®-} 00®, a ground tempera- 
ture of 288®K* 

It has been noted that the ground temperature has very little 
influence, and that a vaKiation of about 20® only slightly 
modifies the corresponding wind profile. In the table can also 
be noted the regions of (byiY) in which the calculation is 
impossible, owing to a sign charge in the denominator enuring 
the iteration. 

9‘. 

The appendix contains a general calculation flowchart and the 
corresponding programs in Fortran language. The programs were 
developed by S. Senotier. 

RESULTS OBTAINED WITH SKYLAB 

t 

The Bureau d 'Etudes de Meteorologie Spatiale of the Meteorologie 
Nationale Frangaise undertoolc an ih-flight experiment on Skylab III, 
which enabled accurate measurements to be taken on a case of 

cloudstreets observed on 19 September 1973 over France. The data 

/ 

f 

obtained through these measurements led to this study. 
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On 19 September 1973, at 1400 T.V. , the astronauts of Sky lab III 
took cumulus alignment pictures over the center of France, over 
a zone covering about 10 square kilometers. The topography and 
general direction of the winds revealed that this is not a case 
of relief waves but truly a case of cloudstreets. The mean direc- 
tion (255®) was taken as the alignment direction, and the 
spacings and Ly were measured from the pictures of SI 90 B. 

The Bordeaux 1200 T.V. radiosounding supplied the values of 
the other parameters required by the model (V^ = 8.2 ms ^ and 240®, 
= 8.7 and 252®, T^ = 295.8, y = 0.0059). 

Figures 1 and 2 represent the pictures of these streets at 
small and large scales. The large-scale pictures were derived 
from SI 90 B, and the small-scale pictures are mosaic montages 
of documents derived from the instrument SI 90 A. 


In the Bordeaux 1200 T.V. radiosounding of 19 September 1973, one 
can note a convective layer thickness of 3500 m (figure 3), 
corresponding to a of 7000 m. is taken as the characteristic 

dimension and the comparison of the six types of cell indicates 
that type No. 3, corresponding to the structure - 7000, 

Ly = 14,000, Ljj = 7000 ra, satisfies the criterion of the mini- 
mization of internal energy dissipation. This type 3 structure 
also' corresponds to the spacings measured from SI 90 B pictures, 
if one adheres to the measurements taken from the largest clouds. 

It can be noted that the inter-street spacing thus corresponds 
to twice the thickness of the convective layer. 
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Type 1, corresponding to a symmetric cell, leads to a numeric 
impossibility. Figure 5, which compares the five calculated 
profiles to the real profile from the Bordeaux radiosounding, 
again shows that the type 3 structure provides the best profile. 
Figure 6 shows the comparison, for each profile, of the total 
<l^iCCi?r»^'^cwthe real profile with the total energy dissipated. 

The error in the calculated profile with respect to the real 
profile is approximately proportional to the energy dissipated, 
confirming the value of the criterion employed. We have thus 
retained the type 3 structure. 

It can finally be stated that, for 19 September 1973, a compari- 
son of the wind profiles supplied by the model and by the 
Bordeaux sounding reveal a relatively slight difference (fig. 7). 
It can also be observed that in both cases, the maximum wind is 
very close to 12 ms \ Moreover, the influence of viscosity can 
account for che reduction in the wind rotation, when going from 
the model to the real wind. 

PROSPECTS 

In conclusion, the result obtained in this particular case, thanks 
to Dkylab, has enabled the construction of a model which provides 

• I 

the wind profile by a sort of t calculation, starting from 

the conditions prevailing at the limits and from a simple morpho- 
logical parameter which is easily measurable. 
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It would be interesting to test this model on a sufficient number 
of radiosoundings by selecting adequate radiosoundings, namely, 
those corresponding to vertical hypothermal structures for which 
it would be possible to determine the thickness of the convective 
layer, and for which wind profile measurements are available at 
the same locations and at the same times, 

.t * 

Should the results of such a comparison prove satisfactory, one 
could then try to apply J:his method to satellite pictures of 
cloudstreets, taking the inter-street spacing as the characteristic 
dimension. The mean thermal gradient can then be derived from the 
surface temperature supplied by VHRR, and the direction of the 
streets themselves could be taken as the wind directions at the 
ground and at the top of the convective layer. 


These verifications only will confirm whether the proposed model 
offers a possibility for operational use. 
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• program RUl 

OIMEMSIOM T 

OIMHMGlOK C 

9I-ME>HSrON-0 

data (DDR (1) 
— nATA(CA5(I) 
REAL LYY 
—REAL- LX,LY, 
COMMON/LIEU 

CtO M M 0 D •/ R-E-SD 

COMHOf(/RARA 
— C-OHHON/SOMD 
DATA PI,G,G, 
—CAS ( 3) =H*SQ 
PAS=1G 

e=-PAS»PAS— 

iPAS=EC 

—READ i- T ( 

1 FORMAT (uAlO 

READ 0- — 

FORMAT (3F1Q 









00 FOO K-1,3 

00 5QG INQ=1, 2 - 

J=J + 1 

J-rEO,~I-)— O-G-T-O-^G-1^ 

300 CONTINUE 

ini l-7si V/f ftcr/i/i 




'wXLf. ~Li/VjHD\(S/ — 

GO TO ?D2 

>0 2-1 F-( R E S { 1 , 1 ) -.-EQ eOE-) -1 = 0 

TAO(MM,IG)=I 
<nn rnfiTTMti'^ 


* 

?00 CONTINUE 

nn n t - 1 1 1 

. 

:>00 IGG (I) = (I-1) *5+5G 


pDT^JT■-7^5 1- - - - . - 


rn format o pas de solution'^, 


S-l+<— b^-^L-V tX=L-2-^,/, 

SIH ,A0X,* 2 - LZ = LY LX = LZ/2-^,/, 

■ — S 1-H— , X-, -=^-^-^1: Z = tY-/^ - L- X- L Z * , / , 

SIH ,40X,^ H - LZ=LY/2 LX=LZ/2^,/, 

SlH~,'4CX r2-8H-5— LZ = IY/(E^^SQRT(2) ) -LX = LZ, V, 

SIH ,4CX,31H 6 - LZ=LY/(2*SaRT(2) ) LX=LZ/2 ) 

P-R-IT-iT— ,11) 

iOl’ FORMAT (IHO, 6 X,1H*, Illiy,/,1H , 117 (IH^) ) 

00- Bb-v—^iM-lr-lS — 

M=1000+100* (MM-1 ) 

3-G0-^RIMT--95h?-t-M-,-(-T-A^-<-MH^-I ) , 1 = 1, 11-) — - 

)00 FORMATdH ,16, 1H + ,11F10) 





■ SUI?ROUTINE MCO£LL(Of ,OIPZl»VEMZl,niR2Z,V(:N?Z,K2,ENERGY, IMPRI) 

— OIMEHSION -UGll-j) iVGtlO) ,UCi(lG» ,V01 (IQ) ,UC2U0» ,VC2(10) ,P(i0»,R(10- 
S),Ul(iC) ,V1(1'3 ) 

— — / P a RA H /T~G — — — 1 L2i*’1Z»0?»GA MM -A A — f 0-irG» G A'MM AtP I-? L-X t b-Y-^-t-Z 

COMMG‘E/RESULT/UC (UOG:) ,VC (40QQ) ,ZZ(A COO) 

REAL LX,LY,LZ»H2,N2,L2 ,LZ? — 

C0MH0N/LIEU/IQATE(4) ,PAStI*^AS 

— c omoH/som/u lt* ( 35-1 > , v r i c 3 e 1 ) , 0 r i ( 351 ) — — 

00 100 J-1,K2 

Z-Z-(-J)~P^-S^(-x)-1r') 

00 CONTINUE 

— G A L L U V ( 0 R r 0 IR Z 1 » V E N 2 1 U R Z 1 , VR Z 1) — 

CALL UV (DR,DIRZZ ,VENZZ,URZZ, VRZZ) 

- — GALOUL - UE- LA-PRE f'lERE-V ALEUR -OU GUE SS- FIELD— 

UG(i) =URZl+PAS*l£-3 

■ — ■' \ AG'(~1~) = VR-ZI*+P-AS“*-1'£'~'4' — ' ~*™ — — ^ ^ 

CALL CALCUKURZlrUGd) ,VRZ1,VG(1) ,K EtP AS»Ult 1) »VK1) »NU) 

Fi NU . €Q . 1 ) - GO- T O-E-dl — — 

PtD-URZZ-Ul (1 ) . 

— Rd-v-vRz-z^i-d > — ■ • — — - 

CALCUL DES AUTRES VALEURS OU GUESS FIELD - UG(N),VG(N) 

BO-EG 0~N^-,iD — ^ ^ ^ — 

UG1 = UG (N) +. UCOl 

VGi=VG(N-) ^ — ■— ^ — ' ^ ^ 

CALL CALCUL(URZl,UGl,VRZl,VGl,K2,PASiUCUN) iVCKK) ,NU) 

OP-l^U G-1- < N)-U-l (• N ) — " ^ — 

DR1=^VC1 (M)-Vl(N) ■ 

UGE^UG4-M-) — 

VGE-VG(N) + . C CO 1 

— — I-F-(-NU r'tQ'i— 1 ) — GO- -T'O — - — 

CALL CALCUL(URZi,UG2,VRZl,VG2,K2,PAS»UC2(H) ,VC2IK) ,NU> 

IR-NU . EO-* 1 )-GG-TG-2-0 1 — — 

DP2=UC2(M)-U1{N) 

BR^ V OE-(-N-)— \MtGN-) — 

UG (M + 1 ) =UG (N) +1 . E-4^ (OR2*P ( N) -DP2»R (N ) )/ (OPl» DR2-DP2»ORl) 

-_„VG < U f i G <iN ) +-1 , E- * ( OP 1 » R ( f4 ) - OR 1 * P- ( K> ) /-( OP dO R 2-DP 2-^OR 1 )— 

CALL CALCUL ( URZ1,UG (N+1) , ORZl, VG (N+ 1) , K2 »PAS i U 1 (N+l ) t Vi (N+1) ,NU) 

I'FT NU^-tiO-i 1) — GO’— T O- ’EO l — ^ — 

P(N+1)=URZZ-U 1 (N-^1) 

RI44+t)-=A/4^32-A/M-^i-^4-) — ^ — 

Wl=P(N-*-l)^^2+R (N + l)^»2 

1 F( W 1 . L T ,-0 -D— G^G-IO-S 0 — — — — — 

’00 CONTINUE 

NU=1 ^ ^ ^ ^ ^ 

VOl EMERGY=1.E06 

GO-I-O ■ 4 . C G ■ ■ — • — - 

50 CONTINUE . ' 

ENERG Y— G~ ^ ^ — — 

L=K2-1 

D0-3G-0-J-2, L — — — — 

DUC= UC(J+1) “2+UC(J)-»-UC(J-l) 

— -gV’G’^ VC-(-U -KI— 2-^^V(H-J )-+ V G ■( J - 1 > — ■ 

... fHERGY=ENERCY*ABS«UCCj.«paC*VClJ.*^^^ OFJOBIAL FA6B IB. / 

21!11N^ I ■ ; OF POOR QOAUrsii - 

IF-( I HPPGl,-EO^)— GO-T^-SO G — — — - — — 

PRINT 8 


*-b"’FOF;MAT <lHi |UCX - NUI^ERICAL A.PPLICAT lOMS-OF LINE A P- THEORY- OF- CLOUDS- 
STREEl S* IMG, GIX,* (WIND pROFT L E ) // ) 

* — PRINT ?1€-,(IDATE(I) ,I = liA) — 

10 FORMATC 1H0,?9X , A A 1 9 i // r 1 H 0 , &AX , ‘D AT A * ) 

PP.i H T— 5 20 L-X , L- Y- » 12 - , V E N Z 1 1 n I R Z 1 , V t N Z Z-, 0 1 R 7 Z , T G , G A H M A 

20 FORMAT ( IH 0 , 4 RX ♦ DP^ 3X^LX * X^L Y*- y X'f-LZ /A gx i F7 , :-:F 1 0 // X + \/ V0*^7 X»DDO» 




-S7X'^VVH*7X^0DH^»JX*TO'f5X»GAMMA^/39X,FR.l,F10,F10, iTFluiFlO.lf F-10 
IF(MU.EQ.l) GO TO 5G1 

-PRINT - 1!', ENERG Y 


ENERGY= '^EIO.SZ) 


11 FORMAT (1H0,//Z ,64X,^M0DEL 

PR-IN T-—1-2 — — — 

12 FORMAT (1H0,E CX ,*HEIGHT» ,7X,* SPEEDS, 6X,+‘DIRECT I ON-^ ,/, IH 


,51X,*(M)-», 


SRXr*' ( M/S > *r 3Xr* < DEG ) ♦ ) 

ERREU1=C 

ERREU 2= E 

00 600 I=1,K2,IPAS 


-fF^SOR-T-(HJG-(-T-) ^^--+--2-+\/G-(-T-) * 2 ) - 


ANG=180 ./PI*ATAN (VC (IJ/UCd) ) 
-00=-0R-ANG— 


ERREUi=ERREUl+SQRT (FF»*2 + VRI (I ) ♦ ^ a-E^FF'^ VRI ( I ) ^CCS ( PI ♦ ( 0 RI ( I ) - DD ) / 


S'180*') ) ' ■■ 

ERHEU2=ERREU2+ AES(FF-VRI (I) ) 

—00^ Amo D-(-Bt}-i-3-6g-i-) 

00 PRINT 1E,ZZ(I) ,FF,OD 

1-3-f ORMAT ( lH-,5-3X', F.5-,-7-X,f 5. 1,&X,F7-) 

PRINT 50E ,£RRE U1 ,ERREU2 

S5-F OP.H A T-{-5-EX-r^IRR'E-y R— !-*■ ♦ E 1 2 . -3 , •/, 3 0-X-, «-eRR£UR~2'^-, E-1-2-.-3 ) 

00 RETURN 

' — PR-I N'T — £0-2 ' 

■02 FORMAT (IHOi^FAS DE SOLUTION*) 

— return — 

END 



S UilROU TIME C AL CU L ( U 1 , U 2 , V 1 1 V 2 r K2 , PA S , U K 2 1 VK2 , M U ) 

COMMON/P APAM/TO ♦ L 2 , M2 , 02 , G A HM A A , 0 , G, G A MM A » PI , LX , L Y , L-Z 

COMMON/RESULT/UC (AQQO) ,VC(4C00 ) ,ZZ(ACOO) 

— — — 

REAL M2,L2 

UC(1)=U1 $ UC( 2) =U2~-S--VC( 1> $-VG(2)=V2 — 

DO 6jG J = 3,K2 

A = G/(TO-GAMMA ^ZZ( J) (M2 + L2) ^ (GAMMA -GAMMAA) 

DENOIJ) =-L2-^UC ( J-1) ^■»2 + M2'»'/C (J-1 ) ^^2 

G-{-2-)-^E^ 0430 

A1 =A^^Q+UC (J“l )/DENO( J) 

A 2 =UG (-J- 1 ) ^4 0 c^-Q- 2 ) +UC ( J - ? ) — 

ai =A^Q^VC( J-i)/(-OEMC(J) ) 

B 2 = V G ( J-- 1 ) ♦ ( 0 2 ♦ G-2 ) 4 V C < J - 2 ) ^ 

UC(J)=A1-A2 

V-G4- J — — ~ 

IF (UENO(J)^DENO(J-l) .LT.D.) GO TO 80 

>-0-0- G OM T I N U E ^ — ^ — ^ 

UK2=UC(K?) $ VK2= VC(K2) 

G-0- TO-^ 0— — — — ^ 

80 NU = 1 

R-trT-U'B'l — - — 

90 MU=C 

R €4 UR-N ^ 



S'UGP.nUTINE UV( DR,DIR, VEN,UR, VR) 

-PI:=3 . 1 A 150 

n=(DR-DTR)*PT/lBO. 

-OP.--V EN^^^OSR-P) 

VR= VCN^SINtO) 

-RETURN ^ 

END 



